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ABSTRACT 

With the aim to investigate the statistical properties and the connection between thermal and non-thermal 
properties of the ICM in galaxy clusters, we have developed a statistical magneto-turbulent model which describes, 
at the same time, the evolution of the thermal and non-thermal emission from galaxy clusters. In particular, starting 
from the cosmological evolution of clusters, we follow cluster mergers, calculate the spectrum of the magnetosonic 
waves generated in the ICM during these mergers, the evolution of relativistic electrons and the resulting synchrotron 
and Inverse Compton spectra. We show that the broad band (radio and hard x-ray) non-thermal spectral properties 
of galaxy clusters can be well accounted for by our model for viable values of the parameters (here we adopt a EdS 
cosmology). 
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I. Introduction 

The most important evidences for non-thermal phe- 
nomena in galaxy clusters comes from the synchrotron 
radio emission diffused on Mpc scales (Radio Halos 
(RH)) observed in a growing number of massive clus- 
ters {e.g., Feretti 2003) and, more recently, from the 
hard X-ray tails (HXR tails) in excess of the ther- 
mal bremmstrahlung spectrum detected in a few cases 
(Fusco-Femiano et al. 2004). RHs are usually found 
in merging clusters (e.g, Schucker et al. 2001) and the 
detection rate of RHs shows a abrupt increase with in- 
creasing the X-ray luminosity and mass of the host clus- 
ters: about 30-35% of the galaxy clusters with X-ray 
luminosity larger than 10^^ erg show diffuse non- 
thermal radio emission (Giovannini & Feretti 2002). 
This observations suggest a strong connection between 
the existence of non-thermal components and the dy- 
namical properties of the thermal ICM. 
A promising possibility to expalain giant RHs is given 
by the presence of relativistic electrons reaccelerated 
by MHD turbulence developed during cluster mergers 
in the ICM (Brunetti 2003; Brunetti 2004 these pro- 
ceedings). Up to now there are no theoretical works 
which investigate if the turbulent reacceleration is able 
to reproduce the statistical properties of RHs in the 
framework of a hierarchical scenario of cluster forma- 
tion. Thus, with the aim to investigate the statisti- 
cal properties and the connection between thermal and 
non-thermal phenomena in galaxy clusters, we have de- 
veloped a statistical magneto-turbulent model. In our 
work, we model the formation of RHs and HXR tails 
in a self-consistent approach which follows, at the same 
time, the evolution of the thermal properties of the ICM 
and the triggering and evolution of the non-thermal 
phenomena in the framework of the magneto-turbulent 
re-acceleration class of models. In particular, we fol- 



low the formation and evolution of clusters of galaxies, 
the generation of merger-driven turbulence in the clus- 
ter volume and the acceleration and time-evolution of 
the relativistic particles, and the related non-thermal 
emission. The goals of our model are: 

1 to check if cluster turbulence generated during 
mergers may be able to drive efficient particle ac- 
celeration processes in the ICM. 

2 to investigate, in the framework of the turbulent- 
acceleration hypotesis, if the hierarchical forma- 
tion process of galaxy clusters can naturally ac- 
count for the observed statistics of RHs. 

In the following we describe our model and investigate 
the point (1); a more detailed discussion is presented 
in Cassano & Brunetti 2004. The second point and 
the Luminosity Function of RHs expected from this 
model are discussed in Cassano, Brunetti, Setti, these 
proceedings. 

II. The magneto-turbulent model 
(a) Cluster formation 

The evolution and formation of galaxy clusters is 
computed making use of a relatively simple semi- 
analytical procedure based on hierarchical Press & 
Schechter 1974 (PS) theory of cluster formation. We 
use the extended PS formalism developed by Lacey 
& Cole (1993) which gives the probability that a par- 
ent cluster of mass Mi at time ti had a progenitor of 
mass M2 M2 + dM2 at some earlier time ^2 , with 
Ml > M2 and h > t2 (e.g., Lacey & Cole 1993, Ran- 
daU, Sarazin & Ricker 2002): 
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Fig. 1. — Example of Merger Trees obtained from 
our Monte Carlo simulations in a EdS cosmology for a 
cluster with a present day mass Mo = 2.5 x 10 ^^M©; 
each cross marks a merger event with M2 > 10^'^ 
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Fig. 2. — Evolution of the thermal energy (solid lines) 
injected in fluid turbulence (dashed lines) in typical 
galaxy clusters. The thin lines are for a cluster with 
present time mass Mo — lO^^M© and the tick lines are 
for a cluster with present time mass Mo = 5 x 10^'' Mq. 
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Where (t(M) is the rms density fluctuation within a 
sphere of mean mass M and bciz) is the critical linear 
overdensity for a region to collapse at redshift z. 
Adopting a suitable change of variable (M with S = 
a^{M) and t with x = Sc{t)) Eq.(l) is replaced by: 



K{AS, Ax)dAS 
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(2) 

Thus the probability that a merger with a given AS 
(i.e. AM) occurs at a given time can be calculated by 
making use of the cumulative distribution of subclus- 
ters masses: 



rAS 

V{< AS, Ax) = / K{AS\ Ax) dAS' = 
"'0 

where erfc() is the complementary error function. 
We employ a Monte Carlo procedure which selects 
a random number, r, in the range 0-1, and determines 



the corresponding value of AS" solving numerically the 
equation: 'P{< AS, Ax) = r. The value of 6*2 of the 
progenitor is given by S2 = S2 + AS and the mass 
by (T^(M2) = 5*2. The mass of the other subcluster is 
AM = Ml — M2 . This procedure is thus iterated until 
the mass of the main cluster drops below a critical value 
or when a maximum redshift of interest is reached. An 
example of merger tree obtained from our procedure is 
illustrated in Fig.l. 

(b) Turbulence in galaxy clusters 

The turbulence in the ICM is supposed to be in- 
jected during cluster mergers. The energetics of the 
turbulence injected in the ICM is calibrated with the 
PdV work done by the infalling subclusters in pass- 
ing through the volume of the most massive one. This 
energetics is given hy Et ~< p >icm vfVt where Vt 
is the volume sweaped by the infalling subhalos which 
is estimated following standard recipes based on Ram 
Pressure Stripping (e.g., Fujita,Takizawa, Sarazin 2003; 
Cassano & Brunetti 2004) and Vi is the impact veloc- 
ity between the two colliding subclusters. Turbulence 
driven by a given merger is assumed to be injected and 
then dissipated within a crossing time, T^ros- We as- 
sume that a fraction rjt of the energy of the turbulence 
developed during these mergers is channeled into fast 
magnetosonic waves (MS waves). Ems = VtEt- 

We use these waves since their damping rate and 
time evolution basically depend on the properties of 
the thermal plasma which are provided by our merger 
trees for each synthetic cluster. As a relevant exam- 
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Fig. 3. — a) a singol merger tree for a cluster with present-day mass of Mo = 1 x lO^^M©, each dot represents 
a merger event with M2 > 10^^ Mq; b) the time-evolution (with redshift) of the electron acceleration coefficient 
{Tacc = Xacc) during the cluster formation; c) the evolution of the electron spectrum at different relevant times (see 
panel) in the same cluster; c) the corresponding radio (synchrotron) and hard-X ray (IC) emission. Calculations 
are performed assuming: Rh = 50Qh^^Kpc, rit = 0.26, T]e = 0.003 and < B >= 0.5^G. 



pie, in Fig. 2 we report the cosmological evolution of the 
thermal energy of galaxy clusters with different masses, 
together with the total energy injected in form of tur- 
bulence in the ICM. The energy in form of turbulence is 
calculated by integrating the contributions from all the 
merger events experienced during the cluster life. The 
thermal energy of the considered clusters, calculated 
assuming the observed M-T relation (e.g., Nevalainen 
et al. 2000), increases from about lO^^erg at z ~ 1 
to a few lO^^erg at the present epoch depending on 
the mass of the cluster. As it should be, we note that 
the energy budget injected in turbulence during cluster 
formation is well below the thermal energy (~ 15% of 
it), which is also in agreement with recent numerical 
simulations (Sunyaev et al. 2003) and with very recent 
observational claims (Schuecker et al. 2004). Finally, 
as reasonably expected, the energy injected in turbu- 
lence calculated with our approach is found to roughly 
scale with the thermal energy of the clusters. 
Under the physical conditions in the ICM, the spec- 
trum of MS waves injected during each merger can be 



estimated as (Cassano & Brunetti 2004): 

/(k) 



rth,e(fc) 



(4) 



where Tth,e is the damping coefficient with ther- 
mal electrons and Ik is the spectrum of the injection 
rate of MS waves (we assume a Kolmogorov spectrum, 
I{k) oc k~^/^ ). We are interested in the energy of 
the turbulence injected in a volume typical of a RH 
[Vh — 47ri?|^/3). The normalization of I{k) is thus 
given by: 



Vh / I{k)dk 

/fen 



(5) 



where /cq is the wavenumber correspondent to the 
maximum scale which is fixed at the stripping radius 
for each merger event (Cassano & Brunetti 2004). 

(c) Particle Acceleration 

We assume the presence of relativistic electrons in 
the ICM which are continuously injected by AGNs, 
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Galactic Winds, and/or merger shocks. Given the cal- 
culated (Sect, b)) spectrum of MS waves and the phys- 
ical conditions in the ICM (Sect, a)), we compute the 
time evolution of relativistic electrons at each time step 
by solving a Fokker-Planck equation including the elec- 
tron acceleration due to MS waves: 



dN{p,t) _ d 
dt dp 
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(dp/dt)c is the Coulomb losses term due to the in- 
teration with the thermal plasma (proportional to the 
density of the ICM) and {dp/dt)rad is the term of ra- 
diative losses due to the synchrotron and IC scattering 
off the CMB photon. For the injection rate of relativis- 
tic electrons we adopt a power law spectrum up to a 
maximum momentum: 



Qe{p) = Kep~ 



(7) 



We normalize the injection rate by assuming that 
the total energy injected in relativistic electrons during 
the cluster life is a fraction, rj^, of the total thermal 
energy of the cluster at z—0. 

The electron diffusion coefficient in the momen- 
tum space due to the interaction with the MS waves 
{e.g., Eilek 1979; Cassano & Brunetti 2004) is given 
by: 



I?pp(p,t) ~ 4.45 7r 



2 «M ^ 



k'W^{t)dk 



(8) 



The acceleration time scale is t„ 



X ^ 2Dpp/p^ 



and thus the systematic energy gain of electrons due to 
MS waves is {dp/dt)acc = XP- The coefficient of elec- 
tron acceleration at redshift z is obtained by combin- 
ing the effect due to the MS waves injected during the 
mergers which contribute to the turbulent spectrum at 
z, one has (Cassano & Brunetti 2004): 



2.23 X 10~^% 
(i?H/500kpc)3 



E 



M1+M2 2.6Mpc\3/2 



.2xl0i5Mo Ri 

(r,/500kpc)2 



(fcr/7keV)i/2 



where Rh is the radius of the RH and Ts is the strip- 
ping radius of the subclumps. 

In Fig. (3) we report an example of our results for 
a cluster with present-day mass Mo — lO^^M©. In 
Fig. (3a) we report the merger history (merger tree) 



which is obtained by using extended P&S formalism. 
Given the merger tree, we calculate the turbulence in- 
jected during cluster mergers and assume that a frac- 
tion rjt is channeled in MS waves {rjt ~ 0.26 is adopted 
in this case) . The coefficient of electron-acceleration is 
calculated considering all the merger events which con- 
tribute to the injection of turbulence at redshift z in 
the considered cluster (Fig. (3b)). The evolution with 
redshift of the electron spectra (Fig. (3c)) is obtained 
solving the Fokker-Planck equation (Eq.JHJl) which ac- 
counts for the relevant energy losses and the acceler- 
(6) ation of the relativistic electrons due to MS waves. 
In Fig. (3) we assume that the energy injected in rel- 
ativistic electrons during the cluster life is a fraction 
{rje — 0.003) of the present energy of the thermal pool 
and calculate the synchrotron and IC emission spec- 
tra at different redshift (Fig. (3d)) from a region of 
Rh ^ 500 kpc h^Q and adopting an average magnetic 
field of < B >~ 0.5 IJ.G. 

The main goal of Fig. (3) is to show that the typical 
observed radio, L_r ~ lO"*" — 10''^ erg {e.g., , Fer- 
etti 2003) and hard X-ray, Lhx ^ lO'*^ - 10^'' erg 5"^ 
{e.g., , Fusco-Femiano et al. 2003) luminosities can be 
obtained in M'^IO^^Mq clusters during major merger 
events, provided that a fraction of the cluster thermal 
energy (of the order of 3 — 5%) is channeled into MS 
waves. The additional requirement is that a population 
of relativistic electrons to be reaccelerated is present in 
the ICM. We find that a total energy injected into these 
electrons of a few 10^^ — 10^"^ times the present energy 
of the thermal pool is sufficient. 
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